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ABSTRACT 

■^j- ' Theoretical models predict that the initial metallicity of the progenitor of a Type la 

(N) | supernova (SN la) affects the peak of the supernova light curve. This can cause a 

C*") ■ deviation from the standard light curve calibration employed when using SNe la as 

OO ' standardizablc distance candles and, if there is a systematic evolution of the metallic- 

, ity of SN la progenitors, could affect the determination of cosmological parameters. 

Here we show that this metallicity effect can be substantially larger than has been 
estimated previously, when the neutronisation in the immediate pre-explosion phase 
. in the CO white dwarf is taken into account, and quantitatively assess the importance 

i' of metallicity evolution for determining cosmological parameters. We show that, in 

principle, a moderate and plausible amount of metallicity evolution could mimic a A- 
dominated, flat Universe in an open, A-free Universe. However, the effect of metallicity 
evolution appears not large enough to explain the high-z SN la data in a flat Universe, 
^ , for which there is strong independent evidence, without a cosmological constant. We 

also estimate the systematic uncertainties introduced by metallicity evolution in a A- 
dominated, flat Universe. We find that metallicity evolution may limit the precision 
, with which f2 m and w can be measured and that it will be difficult to distinguish 

■ evolution of the equation of state of dark energy from metallicity evolution, at least 

from SN la data alone. 

Key words: cosmological parameters - distance scale - supcrnovac: general - super- 
novae: Type la - galaxies: evolution 
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1 INTRODUCTION 

The use of Type la supernovae (SNe la) as standardizable 
cosmological distance candles (Riess et al. 1998; Perlmutter 
et al. 1999; Tonry et al. 2003; Riess et al. 2004; Astier et 
al. 2006) relies on the empirical fact that there is a tight 
correlation between the supernova peak brightness and the 
width of the supernova light curve (Phillips 1993), i.e. the 
fact that, to lowest order, SN la light curves form a one- 
parameter family of curves. The driving parameter that de- 
termines the relation, assuming that the mass of the progen- 
itor white dwarf at explosion is constant, has been shown 
to be the opacity in the ejecta (Khokhlov, Miiller & Hoflich 
1993; Hoflich et al. 1996, Mazzali et al. 2001). This is closely 
related to the quantity of radioactive 56 Ni synthesised in the 
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explosion, which is responsible for the SN luminosity. In re- 
cent years it has become apparent from a larger sample of 
observed supernovae that not all SN la light curves fit into 
this one-parameter family, producing an intrinsic, scatter in 
the Phillips relation (see the discussion in Mazzali & Podsi- 
adlowski [2006] and Benetti et al. [2004]). This immediately 
implies that there must be more than one parameter con- 
trolling SN la light curves. The physical property in the 
progenitor that determines the dominating (first) parame- 
ter still has not been clearly identified. On the other hand, 
from a theoretical point of view it seems unavoidable that 
the metallicity of the original supernova progenitor must 
at least in part be responsible for the intrinsic scatter about 
the mean Phillips relation. Timmes, Brown & Truran (2003) 
showed, using straightforward and uncontroversial nuclear 
physics arguments, that the neutron excess in the immediate 
progenitor white dwarf, which is a direct function of the ini- 
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tial metallicity of the progenitor star, controls the fraction of 
radioactive to non-radioactive Ni produced in the exploding 
white dwarf and hence affects the peak supernova luminos- 
ity, an effect subsequently confirmed by detailed explosion 
calculations (Ropke et al. 2005; Travaglio, Hillebrandt & 
Reinecke 2005). On the other hand, the width of the super- 
nova light curve is ultimately determined by the opacity of 
the ejecta which is mainly a function of the total amount 
of iron-peak elements, i.e. is independent of the ratio of ra- 
dioactive to non-radioactive material (Mazzali et al. 2001). 
Mazzali & Podsiadlowski (2006) demonstrated that this pro- 
duces an intrinsic scatter around the Phillips relation com- 
parable to the observed scatter. An intrinsic scatter in the 
Phillips relation in itself does not necessarily limit the use- 
fulness of SNe la as cosmological probes, as long as the scat- 
ter is non-systematic. However, a physical parameter such 
as metallicity can in principle introduce subtle evolutionary 
effects that could affect the determination of cosmological 
parameters. It is the purpose of this paper to quantify the 
magnitude of such evolutionary effects and to demonstrate 
that these need to be taken into account, in particular when 
trying to constrain higher-order effects such as the cosmo- 
logical evolution of the equation of state of dark energy or 
its variants. 

In Section 2 we show that the effect of metallicity on the 
supernova light curves may actually be significantly larger 
than originally estimated by Timmes et al. (2003), since 
these authors did not consider the additional neutronisa- 
tion in the white dwarf core in the immediate pre-explosion 
phase. In Section 3 we investigate quantitatively the con- 
ditions under which the apparently observed deviation in 
the Hubble relation from a A-free cosmology could be ex- 
plained by evolutionary effects, and in Section 4 we estimate 
the metallicity-dependent corrections when determining the 
equation of state of dark energy in a A-dominated cosmol- 
ogy, as presently favoured by various experiments. Finally in 
Section 5 we discuss these results and possible observational 
tests. 



2 RELATION BETWEEN M PE ak AND Z 

As Timmes et al. (2003) showed, the neutron excess in the 
core of a CO white dwarf at the time of the thermonuclear 
runaway determines the ratio of stable to unstable Ni, which 
in turn determines the supernova peak luminosity. The neu- 
tron excess depends on the abundance of elements that have 
an excess of neutrons, such 56 Fe and 22 Ne. The main source 
of excess neutrons in the interstellar medium from which a 
star is born is 56 Fe which has four extra neutron in each 
nucleus. The neutron excess 1 for 56 Fe therefore is 

7756 = 4X( 56 Fe)/56, (1) 

where we use the notation X(E) to denote the mass frac- 
tion of a given element E. The abundance of 56 Fe is not 

1 The neutron excess, r)i, for a nucleus of species i can be written 
as An Y;, where An is the number of extra neutrons per nucleus 
(i.e. the number of neutrons minus the number of protons) and Y; 
is the number fraction for this species. The total n, summed over 
all species, is directly related to the electron number fraction, Y e , 
according to n = 1 — 2Y e . 



modified until the explosion (any gravitational settling will 
be removed by the growth of the convective core during the 
carbon flash [C-flash]). 

In contrast, 22 Ne is formed predominantly during the 
progenitor's nuclear evolution: during hydrogen burning, the 
CNO cycle converts most of the initial nuclei of 12 C, ie O into 
14 N, which subsequently is converted into 22 Ne by two a 
captures during helium burning. Thus, essentially all of the 
initial abundances of C, N and O are converted into 22 Ne in 
the white-dwarf core, making it the most abundant element 
after 12 C and 16 O 2 . Hence there is a direct relation between 
the 22 Ne contained in the C+O WD and the amount of C, 
N and O nuclei initially present in the star: 

J_ X ( 22 Ne) = ^X( 12 C) + ^X( 14 N) + ^A( lfi O). (2) 

Since 1 O is usually the most abundant element by number 
in star- forming regions, it is the initial 16 O abundance that 
is the most important element for determining the neutron 
excess in the white-dwarf core. 

Timmes et al. (2003) assumed that the neutron ex- 
cess at the time of the explosion was entirely determined 
by the initial abundances of C, N, O and Fe. However, be- 
fore carbon burning runs away and initiates the supernova 
explosion, there is a long phase of low-level carbon burn- 
ing, typically lasting several thousand years, during which 
the excess energy produced by carbon burning is efficiently 
transported away from the burning region by convection, 
delaying the thermonuclear runaway and leading to a grad- 
ually increasing convective core. In this phase, the densities 
are high enough that electron captures are very efficient and 
further neutronise the core material. To model this realisti- 
cally, one needs detailed stellar evolution calculations with 
a detailed nuclear reaction work and a proper treatment of 
the convective Urea process (see Lesaffre, Podsiadlowski & 
Tout 2005). However, with some approximations, we can es- 
timate the effect of additional electron captures onto 22 Ne 
prior to the explosion phase. Each electron capture (or in- 
verse beta decay) increases the number of excess neutrons 
per nucleus by two (since one proton is converted into one 
neutron). On the other hand, at these high densities the in- 
verse reactions (beta decays, positron captures) are strongly 
disfavoured energetically. Therefore, the change in the over- 
all neutronisation of the matter depends only on the total 
number of electron captures/inverse beta decays. 

For this analysis, we consider the reduced network of 
the main nuclear reactions during carbon burning shown in 
table 6 of Arnett &; Thielemann (1985). In this network, 
we further neglect neutron-capture reactions compared to a 
and proton captures. This is reasonable since the main ba- 
sic reactions are 12 C( 12 C, a) 20 Ne and 12 C( 12 C,p) 23 Na which 
produce a particles and protons in approximately equal 
amounts, whereas neutron sources come from much weaker 
reactions. 

We now consider the fate of 22 Ne in this network and 
record the net number of electron captures involved in the 
subsequent chain of reactions. Each electron capture will 



2 Note that the 12 C and 16 O nuclei that exist at the time of the 
explosion have been exclusively produced in the progenitor star 
by the triple-o reaction and a captures during helium burning. 
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indeed add two extra neutrons/nucleus. Neglecting the neu- 
tron capture reaction, 22 Ne has one of two choices: 

• it can either capture an a particle and be locked in 
26 Mg via the chain 22 Ne(a, n) 25 Mg(n, 7) 26 Mg, 

• or it can capture a proton via 22 Ne(p, 7) 23 Na. 

Using analytical rates from the REACLIB (Rauscher & 
Thielemann 2002) data base, we observe that the branch- 
ing ratio of these two reactions is more than 97% in favour 
of the proton capture for all temperatures below 2 x 10 9 K, 
a temperature which is well below the maximum tempera- 
ture that can be expected during the C-flash (Lesaffre et al. 
2006). We can hence safely neglect the alpha-capture reac- 
tion. 

The 23 Na nucleus produced by the proton capture 
can now either capture an electron with a probability q 
to become 23 Ne or capture a proton with a probabil- 
ity 1 — q. In the latter case, the chain of reactions is 
23 Na(p, 7) 24 Mg(p, 7) 25 Al(/3 + ) 25 Mg with one inverse beta de- 
cay. 25 Mg can capture either an electron or a proton, but in 
both cases it leads to one electron capture/inverse beta de- 
cay since 25 Mg(p, 7) 26 Al(/3 + ) 26 Mg. The probability q is not 
easy to estimate since it depends on how the proton-capture 
rate compares to the convective Urea process on 25 Mg. 

A 22 Ne nucleus can hence either capture one electron 
with a probability q or capture two electrons with a proba- 
bility (1 — q). Finally, if we assume that a fraction / of the 
22 Ne is burnt during the C-flash, this leads to an increase of 
the neutron excess of 2 / (2 — q) for each initial 22 Ne nucleus 
which already has two extra neutrons. Summing up these 
contributions, we obtain the total neutron excess before the 
C-flash due to 22 Ne 



?722 = 2[l + /(2-q)]X( 22 Ne)/22. 



(3) 



Based on our present C-flash calculations (Lesaffre et al. 
2006), we estimate that at least 0.02 of the mass of carbon 
is burnt in the C-flash phase, implying that most of the 22 Ne 
may be consumed (i.e. / ~ 1). However, this will generally 
depend on the initial abundances and also the C/O ratio and 
needs to be confirmed with a full reaction network (both / 
and q could in fact depend on the first parameter!). 

Summing up all contributions of metallicity to the neu- 
tronisation from the birth of the star until the start of the 
explosion, we get: 



rjz = V22 + 7?56 



A( 5C Fc) 
56 



+2 [1 + / (2 - q)) 



X( 12 C) X( 14 N) X( m O) 



12 



14 



16 



(4) 



Finally, the neutron excess might change according to 
parameters other than the metallicity, and we write the net 
neutronisation at the time of the explosion as 



n = »?o + vz, 



(5) 



where the subscript denotes the neutronisation for zero 
metallicity. 

Following Timmes et al. (2003), we assume that the 
combustion is fast enough compared to weak interactions 



so that the neutron excess does not change in the process 

rod 

get the mass fraction of unstable 56 Ni 



Allowing only 56 Ni and 58 Ni as products of the burning, we 



A( 56 Ni) 



29ry. 



(6) 



Since the peak luminosity L of a SN la is roughly pro- 
portional to A( 56 Ni) for a given value of the primary pa- 
rameter (e.g. Arnett 1982), we can write 



L = Z X( 56 Ni). 



(7) 



From this we compute the magnitude difference due to 
metallicity effects 



M - Mo 



-2.5 log ! 



29 



29770 



1}Z 



(8) 



where M and Mo are the peak magnitudes for the light 
curves of SNe la of a given primary parameter with (M) 
and without (Mo) metals. 

Setting 770 = (i.e. assuming that the primary param- 
eter does not affect the neutronisation) and adopting the 
solar abundance ratios from Asplund, Grevesse & Sauval 
(2005), we can write equation (6) as 



A( 56 Ni) 



■aZ/Zo. 



(0) 



For Z = 0.02, we obtain a = 0.165 (for q = and / = 1) 
and a = 0.111 (for q = 1 and / = l) 4 . For comparison, with 
/ = (i.e. neglecting the neutronisation during the C-flash), 
we recover a — 0.058, the value obtained by Timmes et al. 
(2003). 

The metallicity dependence of the supernova peak mag- 
nitude relative to a reference model M re f can then be written 

a.s 

M - M rof = -2.51og 10 [(1 - aZ/Z )/(l - aZ Ici /Z Q )] , (10) 

where Z Ic f could, e.g., represent the typical metallicity of the 
SN la sample used in the calibration of the Phillips relation. 
Figure 1 shows this relation for different values of a for zero 
metallicity. Expanding this relation to second order in a, we 
can approximate this relation as 



M - M r , 



1.086a 



Z — Z Y . 



Zo 



1 



1 z 

— a — 

2 



Z 



(11) 



Since the factor in square brackets is not much larger than 1, 
this shows that, to lowest order, the metallicity-dependent 
correction is proportional to both a and (Z — Z Ic f). In the 
following sections, we generally adopt an intermediate value 

3 This will generally not be correct for the central region of the 
white dwarf (~ 0.1 — 0.2Mq) where the main nucleosynthesis 
products are 54 Fe, 56 Fe and 58 Ni (e.g. Nomoto et al. 1984). How- 
ever, this radioactively inert region is not expected to contribute 
significantly to the lightcurve shape (e.g. Mazzali & Podsiadlowski 
2006). 

4 The metallicity normalisation was chosen as Zo = 0.02 since 
this is a typical value adopted in the literature to represent 
"solar metallicity". With the Asplund et al. (2005) composi- 
tion, this implies a logarithmic oxygen abundance by number of 
log(0/H) + 12 = 8.87. In contrast, the corresponding values for 
the Sun are 0.0122 and 8.66 for the Asplund et al. (2005) compo- 
sition mixture, which has a hydrogen mass fraction X = 0.7329, 
and abundance ratios (by number) C/0= 0.537, N/0= 0.132 and 
Fe/0= 0.062. 
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Figure 1. The effect of metallicity on the peak magnitude of 
SNe la, relative to a reference peak magnitude Mo for Z = 0, as 
a function of logarithmic oxygen abundance (by number). On this 
scale, a metallicity of Z = 0.02 corresponds to log(0/H) + 12 = 
8.87 (Asplund ct al. 2005). 

for q of 0.111, but because of this linear relationship, it is 
easy to rescale all the metallicity effects for different values 
of a. 

Hamuy et al. (2000) and Gallagher et al. (2005) have 
previously looked for a correlation between SN la peak mag- 
nitudes and the metallicity of the host galaxy for nearby SNe 
la and found no significant correlation. This indeed confirms 
that metallicity cannot be the dominant first parameter that 
controls the variation of SN la light curves but can only 
be a second parameter. These authors suggested that the 
first parameter is related to age, which in turn could be re- 
lated to the central ignition density (see, e.g., Lesaffre et al. 
2006). Gallagher et al. (2005) also attempted to identify any 
systematic metallicity-dependent deviation from the mean 
Hubble relation and found a weak metallicity trend, but one 
that was not very statistically significant (also see Wang et 
al. 2006). 

Finally, we note that metallicity can modify the prop- 
erties of the exploding white dwarfs in other ways than just 
through the neutron excess; for the single-degenerate model, 
various theoretical studies have predicted or suggested that 
metallicity affects the initial mass range for white dwarfs 
that will ultimately explode (Langer et al. 2000), the C/O 
ratio in the white dwarf (Hoflich, Wheeler & Thielemann 
1998; Umeda et al. 1999) and the accretion efficiency in the 
progenitor phase (Kobayashi et al. 1998; Umeda et al. 1999). 
All of these factors could either enhance or reduce the effect 
of the neutron excess on its own (the possible evolutionary 
consequences of some of these effects have recently also been 
considered by Riess & Livio [2006]). 



3 THE EFFECTS OF METALLICITY ON THE 
DETERMINATION OF COSMOLOGICAL 
PARAMETERS 

As is clear from Figure 1, a systematic variation in metal- 
licity could in principle introduce a systematic effect on the 
observed SN la relation that is comparable in magnitude to 




faint bright 



first parameter (p c ?) 

Figure 2. Schematic diagram illustrating possible evolutionary 
effects when sampling different local and high-z supernovae as- 
suming two underlying parameters. 

the inferred effect that led to the first evidence for an accel- 
erating Universe (Riess et al. 1998; Perlmutter et al. 1999). 
In this section we will therefore ask the question what metal- 
licity evolution is required to mimic a A-dominated Universe 
in a more traditional A-free cosmology. 

The Phillips relation and the first and second 
parameter 

Before investigating metallicity effects, it is important to 
point out that, while the Phillips relation used to correct 
SN la light curves is a one-parameter relation, whether this 
is done using the AM15 method (Phillips et al. 1999) or an 
equivalent method such as the stretch method (Goldhaber 
et al. 2001), this parameter is unlikely to be completely in- 
dependent of metallicity. The calibration of the relation uses 
a local sample of SNe la which includes supernovae over a 
wide range of metallicity. Therefore, the Phillips relation is 
a convolution of whatever the dominant parameter is that 
controls supernova light curves (e.g. central ignition den- 
sity) and the metallicity dependence of the local sample. 
Hence, it is probably not surprising that the local sample 
does not show a clear, systematic metallicity-dependent de- 
viation from the mean Phillips relation (Gallagher et al. 
2005). 

However, this does not imply that there should be no 
metallicity effect at an earlier epoch. As long as there is more 
than one parameter controlling SN la light curves, and as 
long as the supernova sample changes with redshift, system- 
atic deviations from the nearby relation are expected. This is 
illustrated in Figure 2, which schematically shows how sam- 
pling effects can affect the correction function that needs to 
be applied. Note, in particular, that it is possible to have a 
similar distribution of the main observable parameter (e.g. 
AM15) for the local and the high-z sample, but still have a 
systematic bias in the high-z sample 5 . 

5 The consequences of sampling different populations is less 
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Figure 3. Mimicking A-dominated cosmologies. Top panel: the 
evolution of the mean mctallicity (Z) of the supernova sample as a 
function of redshift (z) required to mimic a (O m = 0.3; Oa = 0.7) 
Universe in a A-free Universe with f2 m = 0.0 (empty Universe), 
0.3 (open Universe) and 1 (fiat Universe), as indicated. Q m and 
f^A are the dimensionless densities of matter and "dark energy" , 
respectively, Zo = 0.02, and the present-day Hubble parameter is 
assumed to be 70km/s/Mpc. Bottom panel: the effect of mctal- 
licity on the estimate of the distance modulus for the three cases 
in the top panel, based on equation (10) with a = 0.111 and 
Z re { = Zq. 



In this paper, we will not attempt to correct for an 
intrinsic metallicity dependence in the Phillips relation and 
assume instead that the Phillips relation is just a function 
of the dominant first parameter, independent of metallicity. 
We suspect that this assumption will somewhat exaggerate 
any evolutionary effects due to metallicity. 



Mimicking A-dominated cosmologies 

In Figure 3, we illustrate the metallicity evolution that 
would mimic a flat, A-dominated Universe with f2 m = 0.3 
and Qa = 0.7 in different A-free cosmologies as indicated 
(here fi m is the present matter density as a fraction of the 
critical density and Qa = A/3Hq, where A is the cosmolog- 
ical constant and Ho is the present-day Hubble parameter, 
assumed to be 70km/s/Mpc). These curves were obtained 
by calculating the difference in distance modulus as a func- 
tion of redshift z between the reference A-dominated cos- 



straightforward when the stretch method instead of the AMig 
method is applied, as both local and high-z supernovae are used 
in the calibration (see, e.g., Astier et al. 2006). 
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Figure 4. The mean metallicity evolution of the underlying star- 
forming population as a function of redshift, mimicking a A- 
dominated Universe for the cases shown in Figure 2. The two 
families of curves assume a time delay of 1 and 5 Gyr between 
the star-formation epoch and the supernova, as indicated. 



mology and the respective A-free cosmologies (bottom panel 
of Fig. 3) and then solving equation (10) for Z, where we 
adopted a — 0.111 and assumed that the present-day ref- 
erence metallicity is Zq (i.e. Z Ic f — Zo). For example, if in 
an open Universe with fi m = 0.3, the mean metallicity of 
the supernova sample increased from Zo at z = to about 
2.5 Zo at z — 1, as shown in the top panel, this would exactly 
mimic a A-dominated Universe with Q = 0.3 and S1a = 0.7 6 . 

To relate this metallicity dependence of the supernova 
sample to the metallicity of stellar populations as a func- 
tion of redshift, one also needs to consider that there may 
be a significant time delay between the formation of the su- 
pernova progenitor system and the actual supernova. While 
the typical time delay in one of the most popular single- 
degenerate models is only ~ 1 Gyr, the fact that SNe la also 
occur in old elliptical galaxies suggests that there should also 
be a population of SNe la with relatively long time delays 
(e.g. Cappellaro & Turatto 1988; Branch & van den Bergh 
1993). This could be either due to the tail of the time-delay 
distribution if there is a single dominant progenitor channel 
(e.g. in the double-degenerate channel) or due to a second 
progenitor channel (e.g. Hachisu et al. 1996; for a detailed 
discussion, see Forster et al. 2006). 

In Figure 4, we show the evolution of the metallicity of 
the star-forming population with redshift that produces the 
metallicity dependence in Figure 3, assuming a short time 
delay of 1 Gyr and a long time delay of 5 Gyr 7 . 



An open, A-free Universe with tt m = 0.3? 

As is clear from Figs 3 and 4, in order for the metallic- 
ity evolution to mimic a A-dominated Universe, the metal- 



6 If the reference metallicity were 0.5 Zq, the metallicity would 
have to increase to about 2 Zo , since to lowest order it is change 
in AZ that matters (see eq. 11). 

7 Note that for a given time delay, no SN la would occur beyond 
a particular redshift z max at which the time that has passed since 
the first star formation is equal to the time delay. 
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Table 1. The effects of time delays. 
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Note. — For the given cosmological parameters, to is the look- 
back time at a redshift z = 0.8, z(to + At) gives the redshift 
at lookback times to plus time delays of 1, 3 and 5 Gyr, as 
indicated. 



licity has to increase with redshift, since this makes ob- 
served supernovae intrinsically fainter than is assumed. This 
is at first sight rather counter-intuitive. After all, the global 
metallicity in baryonic matter (stars and gas) has to in- 
crease as the Universe evolves. However, this need not ap- 
ply for the mean metallicity in star-forming regions. In the 
now widely accepted picture of galaxy downsizing (e.g. Treu 
et al. 2005), the most massive galaxies form first and fast 
and have largely finished star formation at a redshift larger 
than 1 (e.g. Cowie 1996; Kauffmann et al. 2003). But the 
most massive galaxies are also the galaxies with the largest 
metallicities (e.g. Kobulnicky & Kewley 2004). Thus, it is 
possible that, as the Universe evolves and star-formation 
predominantly occurs in galaxies of lower mass (with lower 
metallicities), the metallicity in the star- forming regions de- 
creases while the global metallicity increases. Indeed, Pan- 
ter, Heavens & Jimenez (2003) (also see Sheth et al. 2006) 
have claimed that the metallicity in star-forming regions has 
declined for the last 6 Gyr, based on their reconstruction of 
the star-formation and metallicity history of the galaxies in 
the Sloan Digital Sky Survey. However, this method is not 
without its uncertainties (e.g., Mathis, Chariot & Brinch- 
mann 2006). 

More direct evidence that the metallicity in star- 
forming regions may have been higher in the past than at 
the present time comes from the archaeological reconstruc- 
tion of the metallicity history in early-type galaxies (Thomas 
et al. 2005). According to Thomas et al. (2005), the most 
massive early-type galaxies have a metallicity of 2 to 3 times 
solar and formed most of their stars between a redshift of 1.5 
and 2 in low-density environments and a redshift up to 5 in 
the densest environments (see also di Serego Alighieri, Lan- 
zoni & j0rgensen 2006). Thomas et al. (2005) argue that 
the increase in metallicity is mainly due to an increase in 
a elements rather than an increase in iron-peak elements. 
However, since the oxygen abundance is the most impor- 



tant abundance in determining the neutron excess, this is 
directly applicable to our analysis 8 . 

In order to mimic a A-dominated Universe in an open 
fi m = 0.3 Universe, this only requires an increase in the 
metallicity from 1 Zo to ~ 2.5 Zo, well within the maximum 
plausible metallicity range in early-type galaxies. However, 
at a redshift of ~ 0.8 where the SNe la with the largest 
metallicity are required, the most massive galaxies have al- 
ready completed most of their star formation. Therefore un- 
derstanding the time delay between the formation epoch of 
the progenitors of SNe la and the actual explosion becomes 
essential. If there is a systematic time delay of a few Gyr, 
the peak in the metallicity-redshift relation (the solid curves 
in Figs 3 and 4) can easily be shifted into the redshift range 
where massive early-type galaxies form most of their stars 
(also see Table 1). We note that even though the local sam- 
ple of SNe la shows a strong correlation with star formation 
(e.g. Mannucci et al. 2005), implying a relatively short time 
delay of < 1 Gyr (see Forster et al. 2006), this need not 
be the case for the high-z sample where a different popu- 
lation, even one with a longer time delay, could dominate; 
this depends entirely on the details of the star-formation and 
metallicity histories and their effects on the SN la rate. We 
conclude that considering the present uncertainties in the 
metallicity evolution of star-forming regions, a systematic 
metallicity trend as large as required for an fi m = 0.3, open 
Universe cannot be ruled out. 



A closed, A-free Universe with tt m = 1.0? 

In contrast to the fi m = 0.3 case, for a flat, A-free Universe, 
the metallicity of the supernova sample has to increase from 
Zo to about 4Zo at z = 0.8. This is a much larger effect 
than seems plausible, at least for a = 0.111. To produce 
such a large effect, one would probably require the maxi- 
mum value for our estimated range of a (i.e. a ~ 0.165) and 
that the time delays are fine-tuned so as to maximise the 
effect. This would essentially require that all SNe la at a 
redshift around z = 0.8 occurred in the most massive early- 
type galaxies or at least in the bulges of galaxies. This can 
almost certainly already be ruled out from a comparison 
of local and high-z host galaxies (e.g., Strolger et al. 2004; 
Sullivan et al. 2006), which does not show such a dramatic 
trend in host galaxy properties. However, we caution that 
there are still significant uncertainties in our estimate of the 
metallicity evolution effect, due to uncertainties in the value 
of a, the effective metallicity and its calibration to be used 
in high-z galaxies and other cumulative metallicity effects. 
We estimate that this uncertainty is about a factor of 2 in 
either direction. 



8 We note that the underlying evolutionary tracks generally used 
in studies of this type tend to use scaled solar abundances, where 
iron-peak elements contribute significantly to the opacity. The ef- 
fect of using a-enhanced (Fe-deficient) tracks should be to some- 
what increase the inferred metallicity and effective oxygen abun- 
dance, i.e. further enhance the metallicity effect on SN la light 
curves. 
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Figure 5. The effects of mctallicity evolution on the determination of cosmological parameters in a flat, A-dominated Universe. The 
underlying (true) cosmology is assumed to be a flat Universe with f2 m = 0.3, Qa = 0.7 and no evolution of the equation of state (i.e. 
with wq = —1 and w a = 0, where the equation of state is parametrized as w = P/(pc 2 ) = wo + w a (1 — a), and a = 1/(1 + z) is the 
dimcnsionless scale factor of the Universe). The top panels show the mean metallicity evolution in the star- forming component that 
would mimic the indicated deviation from the reference model, while the bottom panels give the effect on the inferred distance modulus 
(the two families of curves in the top panels are for time delays of 1 Gyr [left] and 5 Gyr [right] , respectively) . For each of the three sets 
of models, two of the three cosmological parameters S7 m , u>o and w a are kept fixed, while the third is varied as indicated (SIa = 1 — &m 
in all cases, and Zq = 0.02). 



4 MEASURING THE EQUATION OF STATE 
OF DARK ENERGY 

Independent of the SN la data, there is now ample obser- 
vational evidence for a globally flat Universe from mea- 
surements of anisotropies in the microwave background 
(WMAP; Spergel et al. 2003; 2006) which, when combined 
with other observational constraints (e.g. from galaxy clus- 
tering and weak lensing), provides strong independent evi- 
dence for a flat Universe dominated by dark energy (see the 
discussion and references in Spergel et al. 2006) . It has been 
argued that one of the main challenges in modern cosmol- 
ogy is to understand the nature of this dark energy. Future 
surveys, now under consideration, will attempt to constrain 
the physical nature of this new form of energy by measuring 
its cosmological evolution. It is the purpose of this section to 
quantitatively estimate the systematic uncertainties metal- 
licity evolution is likely to introduce when measuring the 
equation of dark energy in future supernova surveys such as 
SNAP (SuperNova/Acceleration Probe; Aldering et al. 2002; 
also see Riess & Livio 2006 for a complementary study of 
evolutionary effects). 



To describe the time evolution of the equation of state of 
dark energy, we use the parametrization suggested by Linder 
(2003), 

P 

w(a) = — - = w + w a (1 - a), (12) 
pc A 

where a = 1/(1 + z) is the dimensionless scale factor of the 
Universe, the parameter wo gives the present value of the 
equation of state, and w a parametrizes its evolutionary his- 
tory (for a cosmological constant, wo = — 1 and w a = 0). 
For a flat Universe, the Friedmann equation then becomes 

+ (1 - On) a' 3 (i + ^O+^a) e -3u, a (l_a)j _ (13) 

where Ho is the present-day Hubble parameter (again as- 
sumed to be 70 km/s/Mpc). In this parametrization, the evo- 
lution is completely determined by three parameters fi 
and w a (keeping Hq constant). In order to estimate the sys- 
tematic uncertainties due to metallicity evolution, we adopt 
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Figure 6. The evolution of the scale factor as a function of look- 
back time (in units of l/Hg) for a variety of equations of state 
that might be mimicked by metallicity evolution (see Fig. 5). The 
true underlying cosmology is a flat A-dominated Universe with 
f2m = 0.3, wo = —1 and w a = (dashed curve). For the individ- 
ual solid curves, one of the three parameters was varied within a 
range that could plausibly be caused by metallicity evolution (as 
indicated) . 



as a reference model a cosmology with a cosmological con- 
stant with Q m = 0.3, too = — 1 and w a = 0, and then cal- 
culate which metallicity evolution would mimic a deviation 
from this reference model. Specifically, we keep two of three 
cosmological parameters fixed and vary the third. The re- 
sults of this exercise are presented in Figure 5, which in the 
top panels shows the necessary metallicity evolution in the 
star-forming material (for supernova time delays of 1 and 
5 Gyr) as a function of redshift and in the bottom panels 
the implied effect on the estimated distance modulus. 

In order to estimate the systematic errors introduced 
by metallicity evolution, one in principle needs an estimate 
for the range of metallicity variation that is "reasonable". 
Following the discussions in the previous section, we make 
the simple and conservative assumption that a variation of 
Z from a few tenth Zo to ~ 2 Zo could be plausible. We can 
then directly read off rough estimates for the systematic 
errors for each parameter from Figure 5. We obtain Af2 m = 
-0.1/ + 0.05, Aw = -0.4/ + 0.2, and Aw a = -2/ + 1. With 
some prior assumptions about the metallicity evolution, one 
would obviously be able to constrain these errors further. 

To illustrate the uncertainty this could introduce in de- 
termining the nature of dark energy through its time evo- 
lution, we plot in Figure 6 the scale factor as a function of 
lookback time (similar to Fig. 1 in Linder 2003). Here the 
dashed curve represents the evolution of the scale factor for 
the assumed underlying cosmology (a flat Universe with a 
cosmological constant), while the solid curves show the evo- 
lution of a for cosmologies that could plausibly be mimicked 
by metallicity evolution. Comparing this figure to Figure 1 



in Linder (2003), we note that there is a large overlap of 
these curves with the curves produced by different physical 
models of dark energy. This implies that it will be difficult 
to distinguish these different models, unless metallicity evo- 
lution effects are taken into account and are corrected for. 

Our main conclusion therefore is that the possible sys- 
tematic uncertainties in fi m and wo caused by metallicity 
evolution could be quite significant, but may be manageable 
(in particular, if further constraints on the metallicity evo- 
lution can be employed). In contrast, the uncertainty in w a 
seems uncomfortably large, implying that it may be very dif- 
ficult to differentiate time evolution in the equation of state 
from time evolution in Z using SN la data alone. 



5 DISCUSSION AND CONCLUSIONS 

Our main conclusion is that metallicity should affect the 
SN la calibration method; this relies only on reasonably 
straightforward and well understood nuclear physics. As 
long as there is any significant cosmological evolution in the 
metallicity, which a priori is to be expected, this has to af- 
fect the determination of cosmological parameters at some 
level. The main remaining uncertainty is the level at which 
this becomes important. 

In our estimates of systematic errors, we adopted a — 
0.111, roughly twice the estimate of Timmes et al. (2003), 
though a could possibly as large as 0.165. It still remains to 
be shown theoretically that the metallicity effect is as large 
as these estimates suggest. This will require detailed stellar 
evolution calculations of the pre-explosion phase that model 
the convective Urea process realistically (Lesaffre et al. 2005 
and references therein) combined with up-to-date explosion 
modelling (e.g. Ropke et al. 2005), work that is presently in 
progress (F. Forster, et al.). In general, since the magnitude 
corrections scale linearly with a, at least to lowest order (see 
eq. 11), all estimates of systematic errors in this paper can 
be rescaled with whatever the appropriate value for a will 
turn out to be. 

It would be particularly useful to be able to constrain 
the value of a observationally. This is quite a challenging 
task even for the local SN la sample, since it is clear that 
metallicity cannot be the primary parameter controlling SN 
la light curves (Hamuy et al. 2000; Gallagher et al. 2005) 
and without a better understanding of the dominant pa- 
rameter and how it affects the basic light curve correction 
method, it will be difficult to isolate the metallicity contri- 
bution. On the other hand, since the predicted metallicity- 
dependent correction factor is, to lowest order, a linear func- 
tion of metallicity, with a sufficiently large sample of SNe la 
with well determined host galaxy metallicities, one could in- 
clude a linear metallicity term in the Phillips relation and 
attempt a two-parameter calibration. This could in princi- 
ple allow a measurement of the magnitude of the metallicity 
effect (i.e. constrain a) and could then be used to correct 
for metallicity effects at high z 9 . 

9 Obtaining the metallicity of the SN la progenitor directly is 
probably not feasible, and therefore one has to rely on a statis- 
tical measure of the metallicity using the host galaxy metallicity. 
This itself is not without problems, considering, e.g., that there 
may be a substantial time delay between the formation of the SN 
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More progress is likely to come from a systematic com- 
parison of the properties of nearby and high-z supernovae, 
one of the major objectives of both the ESSENCE project 10 
and the SuperNova Legacy Survey (SNLS) 11 . First results 
from the SNLS (Astier et al. 2006; Sullivan et al. 2006) in- 
dicate that the nearby and high-z samples are similar, but 
that there may be some differences in properties. While at 
present these are only of marginal statistical significance, 
this may change as the sample increases. 

The evolution of galaxy properties, in particular the 
metallicity in the star-forming regions, clearly plays an im- 
portant role in the theoretically predicted evolution of SN 
la properties. There is a general consensus that there has 
been a significant evolution of galaxy properties since a 
redshift of 2-3: star formation first occurred in the most 
massive galaxies and since then has moved increasingly to- 
wards lower-mass galaxies. Observationally, this is reflected 
in star formation predominantly occurring in ultraluminous 
infrared galaxies (ULIRGs) at a redshift z > 2, in luminous 
infrared galaxies (LIRGs) at a redshift z ~ 1 and starburst 
galaxies in the local Universe (see, e.g., Fig 10 in Pcrcz- 
Gonzalez et al. 2005). If this galaxy trend is also correlated 
with a trend in metallicity, it would exactly produce the 
type of evolutionary effect that could affect the determina- 
tion of cosmological parameters. A better understanding of 
the metallicities in these different types of galaxies is there- 
fore clearly needed. 

In summary, metallicity effects should be expected and, 
at some level, will affect the determination of cosmologi- 
cal parameters. Uncertainties in the metallicity evolution in 
star-forming galaxies and in the intrinsic properties of SN la 
progenitors (in particular, their time delays) will limit the 
precision with which the equation of state can be measured; 
it will be particularly difficult to distinguish an evolution 
of the equation of state from evolution in SN la proper- 
ties. As the statistical errors in SN la surveys and comple- 
mentary surveys (CMB, galaxy clustering, weak lensing) de- 
crease, this may first show up in statistical discrepancies in 
the determination of cosmological parameters using different 
methodologies. 

Our estimates show that metallicity evolution could 
possibly be large enough to mimic a A-dominated Universe 
in an open Universe, but is unlikely to be large enough to 
mimic a preferred flat Universe. Nevertheless, considering 
the fundamental shift in our view of the physical world adop- 
tion of a cosmology dominated by dark energy would bring 
about, we emphasize the importance of taking such system- 
atic effects into account. One can hope that with a better 
understanding of such effects, one will be able to correct for 
these and ultimately obtain a reliable calibration of cosmo- 



progcnitor and the explosion and that different metallicity meth- 
ods have to be employed for different types of galaxies, all with 
their own calibration issues. In particular, it will generally be dif- 
ficult to obtain a reliable metallicity estimate for the progenitors 
of SNe la in early- type galaxies, since the formation of the super- 
nova progenitor may be associated with a relatively recent merger 
event, in which case its metallicity need not bear any relation to 
the metallicity of the dominant old population (e.g. Nolan ct al. 
2006). 

10 See http://www.ctio.noao.edu/~wsne/indcx.html. 

11 http://cfht.hawaii.edu/SNLS. 



logical parameters, where systematic uncertainties are min- 
imal. 
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